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Abstract 

A numerical method is presented to determine electromagnetic shielding 
effectiveness of rectangular enclosure with apertures on its wall used for input and output 
connections, control panels, visual-access windows, ventilation panels, etc. Expressing 
electromagnetic fields in terms of cavity Green’s function inside the enclosure and the free 
space Green’s function outside the enclosure, integral equations with aperture tangential 
electric fields as unknown variables are obtained by enforcing the continuity of tangential 
electric and magnetic fields across the apertures. Using the Method of Moments, the integral 
equations are solved for unknown aperture fields. From these aperture fields, the 
electromagnetic fields inside a rectangular enclosure due to external electromagnetic sources 
are determined. 

Numerical results on electric field shielding of a rectangular cavity with a 
thin rectangular slot obtained using the present method are compared with the results obtained 
using simple the transmission line technique for code validation. The present technique is 
applied to determine field penetration inside a Boeing -757 by approximating its passenger 
cabin as a rectangular cavity filled with a homogeneous medium and its passenger windows 
by rectangular apertures. Preliminary results for, two windows; one on each side of fuselage 
were considered. Numerical results on electromagnetic field penetration through passenger 
windows of a Boeing -757 at frequencies 26 MHz, 171-175 MHz, and 428-432 MHz are 
presented. 
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1.0 Introduction 

Apertures of various sizes and shapes on a metallic enclosure are used for reasons such as 
input and output connections, control panels, visual-access windows, ventilation panels, etc. 

Since these apertures at appropriate electromagnetic (EM) frequencies behave as very efficient 
antennas, they also become sources of electromagnetic interference (EMI) problems for both EM 
emission and susceptibility. It is important to know the EM shielding effectiveness of these 
enclosures in presence of the apertures. The EM shielding effectiveness study may also help in 
locating these apertures at proper places to reduce the EM emission or improving the immunity 
of electronic components present inside the metallic enclosnre[l]-[5] 

Shielding effectiveness can be calculated using numerical or analytical methods. Numeri- 
cal methods such as the Finite Difference Time Domain (FDTD) method (6], and the hybrid 
Finite Element/ Method of Moments [7] can model complex structures inside enclosure but often 
requires large computing time and memory. For electrically large size enclosures and apertures 
these methods , though more accurate, become difficult for designers to use to investigate effects 
of EM shielding on the design parameters. 

Pure analytical formulations described in [3] -[5] even though provide a much faster means 
of calculating shielding effectiveness are based on various simplifying assumptions whose valid- 
ity may be questionable at high frequencies. 

In the present paper a method which is suitable for large but regular shaped enclosure and 
apertures is described. Using the equivalence principle as described in earlier papers [1], [2], the 
apertures are replaced by equivalent magnetic current sources. By matching the tangential elec- 
tromagnetic fields across the apertures, coupled integro-differential equations with the magnetic 
currents at the apertures as an unknown variables are obtained. The coupled integro-differential 
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equation in conjunction with the method of moment are then solved for unknown magnetic cur- 
rent amplitudes. One of the advantages of the present approach is that it gives an appropriate 
model for EM shielding effectiveness studies and avoids some of approximations or assumptions 
which are made in analytical formulations to arrive at closed form formulas for EM shielding. 
Furthermore, by representing the magnetic current over the apertures in terms of entire domain 
basis functions, three to five number of unknowns per aperture are required for convergence. The 
present method, therefore, requires considerable less computer time compared to other numerical 
methods. The CPU time taken to determine EM shielding for a rectangular cavity illuminated 
through two rectangular apertures with three entire domain modes on each is around 23.00 sec- 
onds on a SGI Origin 2000 machine. 

The remainder of this report is organized as follows. A general formulation of the prob- 
lem of a rectangular cavity with a series of rectangular apertures on its broad walls which are 
illuminated by plane wave at normal incidence is presented in section 2. Section 2 describes: (1) 
use of cavity Green’s function to determine electromagnetic fields scattered due to the apertures 
inside the cavity, (2) use of free space Green’s function in plane wave spectrum form to determine 
scattered field outside the cavity due to apertures, (3) setting of integral equation and use of 
method of moments to determine matrix equation with aperture fields as unknown variables. 
Numerical results on electric field shielding due to horizontally and verically polarized incident 
waves for various sizes of rectangular cavity and rectangular apertures are given in section 3 along 
with earlier published data for comparison and validation. Tire report concludes in section 4 with 
remarks on the advantages and limitations of the present method. Section 4 also gives a recom- 
mendation for future work to be persued. 
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2.0 Formulation of Problem 


Figure 1 shows a geometry of rectangular enclosure with rectangular apertures on its two 
opposite walls. These apertures may be used for input and output connections, control panels, 
visual access windows, etc. Since the rectangular apertures are very efficient electromagnetic 
radiators, they also act as sources of electromagnetic interference both for emission and suscepti- 
bility. The electromagnetic energy radiated by conducting wires and loops inside the enclosure 
may escape through these apertures and cause an interference with another electronic 


a 

◄ ► 



Figure 1 Geometry of rectangular enclosure with rectangular apertures 
illuminated by a plane wave at normal incidence. 


components. The electromagnetic energy from sources outside the enclosure will get through 
these apertures into the enclosure and may cause interference with electronic circuits/components 
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present inside. In this report the latter problem is considered. To determine electromagnetic 
shielding effectiveness of the enclosure shown in figure 1 due to external EM radiation, a plane 
wave at normal incidence is assumed. This assumption is valid for worst case study. The tangen- 
tial electric field induced on the apertures due to the incident field may be represented by 


p. q 




XX,XX^ P ,cos (fffh 
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W r ( 2 
+ x- jc^Y) sin 
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where (x cr , y cr ) are the coordinates of center of r aperture, U rpq and V rpq are the 
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th th th 

unknown voltages of pq mode on r aperture, W r , L r are the width and length of r aperture, 
respectively, x, y are the unit vectors in x, y directions. In equation (1) it is assumed that there 


are R apertures on the z = 0 plane. The unknown modal amplitudes U * 0 and V rpq * Ofor 
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and A rnn , B rnn are the unknown modal amplitudes corresponding to the apertures on z = c 

L r L r 

plane and they satisfiy A rpq ^ 0 , and B rpq ^ 0 for x cr - — < x < x cr + — , 

W W 

y cr — — < y < y cr + — , and A rpq = B rpq = 0 otherwise. The problem of EM coupling to a 

rectangular enclosure through apertures therefore can be split into internal and external regions as 
shown in Figure 2. Assuming the z = 0 and z = c planes are infinite, the entire problem can be 
split into three regions: region I for z < 0 , region II for 0 < z < c , and region III for z > c . It 
should be noted that because of the presence of the infinite ground planes at z = 0 and z = c , 
the apertures on the z = 0 plane do not couple to the apertures on the z = c externally. 

The internal problem consists of a rectangular volume enclosed by the enclosure and illu- 
minated by various equivalent magnetic current sources. The external problem consists of mag- 
netic current sources backed by infinite ground plane for each wall of the rectangular enclosure. 
The unknown amplitudes V , U r , A , and B appearing in (2) and (3) are determined by 

setting up coupled integral equations. The expressions for the EM fields required to form integral 
equation are determined in the following sections. 
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Wave 


Figure 2 Geometry of equivalent problem with apertures replaced by 
magnetic current sources. 


2.1 Electromagnetic Field Outside Enclosure 

The total electromagnetic field outside the enclosure is obtained by superposition of the 
field due to incident wave and the scattered field due to apertures. In the following sections the 
EM field due to incident wave and the scattered field due to aperture are determined. 

2.1.1 Electromagnetic Field Due to Incident Wave 

The incident field with a time variation of e ]m may be written as 

Hi = kl ' f ' 

= (0,|//,| cos(a 0 ) + <j>,;|//, : | sin(a 0 ))e ^ (4) 
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where ki • r = -£ o sin(0 ( )(;rcos((|) ( ) + ^sin((j) / )) -z£ o cos(0,) , k 0 is tlie free-space wave 
number, and (0,-, 4>,) are the angle of incident plane wave. From equation (4), the x-, -y, and z- 
components of the incident magnetic field may be written, respectively as 
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For nomial incidence, with a 0 = 0,0, = 0 , and (]>, = 0 , the incident field in the z = 0 is 


given by H xj = \Hi\, H yj = 0 , and H zj = 0 . 

2.1.2 Electromagnetic Field Scattered due to Apertures: 
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Consider the r aperture on the z = 0 plane, the EM field scattered due to M,-\ can be 
obtained from [8] 
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In expressions (12)-(17) § rpqy is the Fourier transform of O and \|/ is the Fourier trans- 
form of 'P . The expressions for <|> and \|/ are given in Appendix III. 


th — ) 

Consider the r aperture on the z = c plane, the EM field scattered due to Mr 2 can be obtained 
following the same procedure 


R 


hi 

X - 

r = 1 P-q 


= V w- .Ik.dk 

4k 


~jk z \(z-z)\ jk x + jk y 

p ■ / 

rpqy \r' 


(18) 


— oo — oo 


B. 


K - II F 


-Azi(z - z’)i + vF? 

dk x dk y 


r=\pq 


— oo — oo 


(19) 


E III = ^ ^ 1 J J ( A rpAp</> + Kp^frpq^yl e jk x x + ^ 

= : pq 4K 


X y 


— OO — oo 


(20) 


H UI = I I ~ G)e ° Brp(l 

r = 1 Pq 4tT£o 


oo oo 2 2 

f f ~jk 2 \(z-z')\ ( k 0~ k x) jk x x + jk y 

e \ t e dk„dk 


rpqx £ 


* y 
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+ xx^% 

= i pq 4k ko 


CO oo 

j j 


n-pqy £ 


* y 


( 21 ) 


H 


m 


OO OO 2 2 

COEn^,-,,., r r -jk 2 \(z - z')| A (^0 ~ k y) Jk X x + jk y y ^ 

^ /jii f ^ 


- XX^¥'JJ 

/- = i pq 4k ko _ 


rrpqy £ 


* \y 


R 

+ XX 


7‘^eV/ j j Gill/-** ik ’ y dk.dk.. 


2 ^ 

= i pq 4n~ko 


pqx k 


x y 


(22) 


MI 


D oo oo 

-C08, 


XX /•// 

= 1 pq 4 71 &0-OO-C* 


-y>J(z-z')| , . jk X x + jky 

e ( B rpqVrpqx k x- A rpq < $>rpqy k y) e dk x dk 


y' 


2.2 Electromagnetic Field Inside Enclosure 

The equivalent magnetic currents, present on the apertures of an enclosure radiate electro- 
magnetic fields inside the enclosure. The total EM field at any point inside is obtained by a super- 


r/z 


position of fields due to each equivalent magnetic current source. Consider the r aperture, the 


th 


electromagnetic field inside the enclosure due to the r aperture is obtained from 


->ii i ->n 
E = — -VxF 
e o 


(23) 


— /CO 9 ->// ->// 

H = -4f(*5F + VV*F ) 

b 1 

k 0 


(24) 


where the assumed variation e 7tor has been suppressed. The electric vector potential appearing in 
(23) and (24) satifies the inhomogeneous wave equation 


V-F (x, y, z ) + k~ 0 F = -e 0 M y, z ) 


(25) 


If G m (x, y, z/ x\ v\ z ') is the dyadic Green’s function for the rectangular enclosure for a unit 

dyad I(x.,y\z) = xx. + yy + zz. inside the enclosure, then the electric vector potential F can 
be written in the form 
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F (x,y,z) = JJ J G rn (x,y,z/x',y',z')»Mr(x',y',z’)dx'dy'dz' 


Source 


Substitution of (26) in (25) yields 
V 2 G m + k 2 () G rn = -E 0 Ib(x. - x')b(y - y')8(z - z') 


Since the apertures are in the xy plane and M r can have both the x- and y-directions, 
I(x', y\ z') = xx + yy , hence equation (27) can be written in component forms as 


v2G mxx + k o G m XX = -e 0 8(x-x')8(y-y , )8(z-z') 


V 2G myy + k o G m yy = -e 0 8(x- x')b(y - y')b(z - z’) 


2.2.1 Electromagnetic Fields Scattered Due to Apertures on z = 0 Plane 


Considering the x-component of the magentic current and using the proper boundary con- 


ditions, the solution of equation (28) can be written as 


„ £ o s om £ o« . ( mnx'\ f nKy'\ C 0 S (^i( z ~ c )) . ( mnx\ (mty\c.. 

G rr = > 7 ; — sin cos — - — - — sin cos — 5(z ) 

** ^ kj ab \ a ) \ b ) sin (k,c) \ a J \ b ) 


for the aperture located in the z = 0 plane. In expressions (30) !<•>= I l<.) 


m = On = 0 


, |,2 ( mu '] 2 ( nn \ 2 r (rm \ 2 . , . ( mnS 2 ( mi \ 2 ,2 _ 

*' = -[t + (t and *' = t + T -*° for 


k~Q < ’ £ om = 1 form = 0 and e 0w = 2 for m ^ 0 . Substituting (30) into (26) 


->11x0 -> 

the electric vector potential F due to the x-component of M r on the plane z 


= 0 is obtained 
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oo 


JIx 0 


8 0 e 0m 8 0n 


X “U “um“un / / / . . • ( ( nny\ 

. , , — -cos(£,(z-c))sin cos — ^ 

k , absm(kfC) 1 \ a ) \ b ) 


y, 0 ) • 1 ) sin cos ^^-Jjdx'dy' 


( 31 ) 


-,11x0 


The total electric vector potential due to all the apertures is then obtained by superposition as 

;- = IP. <7 m , n 

-> 


<7bsin(^ / c) 7 v tf ) V b 

. ( mnx'\ ( nnV'' 

~'"l Icos 1 


ff( v t', -£<?(*', y') • x)sin( - — — jcosf / 3-3- ] dx'dy' 

\ a ) \ b ) 


(32) 




Since 'fVp<y(.r', /) = Jtr'F ^(jc', /) , expression (32) can be written in component form as 

R 00 ~ ~ ~ 

. (mnx\ (nny' 


IIxO v-i v 1 ~uw~u» , .. . I run a \ / mi 

F r = > > U rm > . — — -cos(£,(z-c))sin cos - 7 ^ 

* ^ ^ jfc absm(^c) v /v \ a ) \ b 

- - 1 n " m,n 


8 0 8 07/7 8 0» 


;-= IP- <7 


(33) 


The total magnetic field inside the enclosure is then obtained from (24) and (33) as 


£ 0 8 0m 8 0« 


k~ ‘ 

1 - 


0 r = 1P’<? 


m , n 


2 ( mn\ 2 \ . ( mnx\ (nny 

I7absin(A- / c)r o_ lvJ M~J C0S l~r 


^f X0 = X ^rpq'L-r^ ~u71,7 ^ k Q-\ — I l Sin l — l C0S l 


cos(^ / (z-c))/ ; . 


pqrnnx 


(34) 


// 


//xo - 7 C 0 ^ ^ 77 _8 o 8 077? 8 0 t7 7 U 7 tf 7 t 7 t^ ( mnx\ . (nny 

— X X^Xi7 gfcsin(t;C) vbyJ c ° s l— H it 


c 0 7 - = IP. <7 


m, 77 


cos(^ / (z-c))7 


r pqrnnx 


(35) 


„/7x0 _ - 7 O) V V 77 V £ ° 

2 ,2 U rpq Zj k. 


e 0m e 0n mn . 


c 0 7 - = IP. <7 


... f mnx's ( nny 

1 u • n \ (-*/) cos cos \-f- 

kj aosin(^ / c) a \ a ) \ b 


sin(£ z (z - c))I 


r pqrnnx 


(36) 
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In (34)-(36), I rpqmnx is equal to JJ v F rp ^(^ t , y)sin^-^^- jcos ^^jdx'dy' and can be 

<7 

expressed in closed form (see Appendix I ). 

Likewise, considering the y-component of the magnetic current and using the proper 
boundary conditions, the solution of equation (29) can be written as 


£ o 8 Om £ o« fmjuA . (nny'\ cos (ki( z ~ c )) (mnx\ . (nny 

n I I nm I I I I nui I *■' 


rnyy 


I7^ c0s [^j sm [-rj sin(v) cos [^J sin [^J 5(z) (37) 


for the aperture located in the z = 0 plane Substituting (37) into (26) the electric vector poten- 

->IIyO — > 

tial F ’ due to the y-component of Mr due to apertures in the z = 0 plane is obtained as 


pllyO _ 'ST' £ 0 £ 0m £ 0« 

t y 2-4 


trn yjn // / \\ ( . f HIZV 

1 u • n — tCos(^(z-c))cos sin -rf- 

k l oosm(k / c) 1 \ a J \ b 


JJ[(M,-(y, y') • y)cos(^^^jsm^-^jjrfyrfy’ 


(38) 


The total electric vector potential due to all apertures on z = 0 is then obtained by superposition 
as 

R 

My o 


£ 0 £ 0w £ 0n 


= y y v y 

y rpq ^ k, absin(k r c) 

r&\lP><l m,n 1 ‘ 


n , (mnx) . (nny 

cos(A: / (z-c))cosl — — Isinl 


jj(<Z>rpq(x', y') • 3 > )cosf^-lsinf^ )dx'dy' 


(39) 


Since Q?rqp(x', y') = —y<S> r/tn Jx',y') , expression (39) can be written as 

’ Hr J 

R 


My 0 


= X I-%I 


£ 0 £ 0w £ 0« 


;• = IP- <7 


m, n 


// / w . fwry' 

7 • ,, — :COs(k / (z- c))cos sin —■ 

absmikjC) 1 \ a J \ b 


|J« > rp W U'.3’') co s(^ £ )sin(^j*'a!y 


(40) 


The total magnetic field inside the enclosure is then obtained from (24) and (40) as 
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CO 


X X-’V/X 


£ 0 £ 0m £ 0« 


‘•O r = IP’ c l 


k[ absm^kjc)^ a ) b 


mn\nn . (mnx\ nny 
—sin cos — 


cos (kjiz-c))! 


rpqrnny 


XX-^X-r 


£ o £ Ow £ o» ( ,2 (nny] ( mnx ) . ( nny 
kj absm(k { c)\ 0 v b J yv a J sm [ ^ 


c 0 r = IP’ 1 


cos (kj(z-c))I 


£ o £ Om £ o» nn, 


uj T7 U Urn \)n nil, 7 . 

2 “ X Yj~ V rpq Yjy~ ab sin(^c) ~b z) C0S 


c 0 r = 1 P’4 


mnx\ nny 

COS — r^- 

a J \ b 


sin (k[(z - c))I r 


rpqrnny 


where I rpqmay appearing in (41)-(43) is equal to jj®rpqy(*> >’’ ) c o s J J sin J dx dy 

q 

and can be expressed in closed form (see Appendix II). 

2.2.2 Electromagnetic Fields Scattered Due to Apertures on z = c Plane 

Considering the x-component of the magnetic current and using the proper boundary con- 
ditions, the solution of equation (28) can be written as 

m, n 

for the aperture located in the z = c plane. Substituting (44) into (26) the electric vector poten- 

-yllxc — ^ 

tial F due to the x-component of Mr on the plane z = c is obtained as 

oo 

r iixc ^ _£ o E 0m E 0n ,, , . ( mnx\ (nny\ 

F r = > t— t— 7 — cos(£rz)sin cos ~r~ 

x ^ k[ absm(kjc) 1 \ a ) \ b J 


JJf (M r (x', y, c ) • x.) sin f ml ^ X yoJ n%y 


dx'dy' 


The total electric vector potential due to all apertures is then obtained by superposition as 


£ 0 E 0m E 0n 


y^nxc n y-' w u w Um w U« \ • ( WIKX\ ( /271V 

F x = > > B rna > — 7 — -cos(krz)sm\ cos — - 

" " p q *- i k, absm(k,c) 1 \ a J \ b 


r = IP’ f l m , n 


kj absinikjc) 
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[ \C¥ r pq(x', y\ c) • 1) sin j cos 7 2_A. dx'dy' 

JJ \ a J \ b J 


Since ' ¥,-pq(x y\ c ) = — jr'P (jt', y') , expression (46) can be written in component form 


F IIxc = YY (_B ) Y g ° £ ° w£ °” cosf/c z)sinf m7tX lcosf n7l}; 

* " " rqp ^ k, absin(k,c ) ^ v a J v b 

r = IP’ 1 m,n 1 7 


( mnx \ (nny 

Vrpqx^^y)* m[— Jcos' 


dx'dy' 


The total magnetic field inside the enclosure is then obtained from (24) and (47) as 


jllxc _ —J CO 


X ’t-' , n ■. ^0 £ 0m £ 0» ( .2 ( HI K \ \ . ( fHKX\ ( nny 

1 17 afcsin( ^ c) (^0-(- ) M — M — 


"0 r = IP’ 4 


— cos It 


cos (kjz)I 


■pqrnnx 


P Q r = \P> m , n 


_£ o £ 0 m £ 0 n mn( nn\ ( mnx\ . (nny 

kj absin(kjc) a riT J C ° S lV lib 


cos(k,z)i 


r pqrnnx 


zln y y (~b t y £ ° £ ° m&o ” mK ( k i cos f 

,2 2-i 2-i rqp 2-i f C{ absinikjc) a 1 1 n J i h 


c 0 r = IP’ 4 


sin(^ / z)7 


r pqrnnx 


Considering the y-component of the magnetic current and using the proper boundary con- 
ditions, the solution of equation (29) can be written as 


r = Y 

myy Lu h 


_£ 0 £ 0 m £ 0 n (mnx'} . (nny'} cos (ki z ) ( mnx} . (nny}?, , 

j;—‘ cos l— j sm l^Jii^ cos l— J sin lirJ 8(z - c) 


for the aperture located in the z = c plane. Substituting (51) into (26) the electric vector poten- 

->IIyc — ^ 

tial F due to the y-component of M r due to apertures in the z - c plane is obtained as 


j^iiyc ^ £ o £ Om £ o« , mnx\ . (nny 

F x = > - . — -cos(k f z)cos sm \ 

y ^ kj absm(kjc) 1 \ a J Kb 
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JJ 


(M r ( x', y\ c) • y)cos(^p- Jsinp^. Jjdx'dy' 


(52) 


The total electric vector potential due to all apertures on z = c is then obtained by superposition 
as 

R 


u yc = y y A V — — — — cost k zlcosf — 1 sin f ^ 1 
y 2-i 2-i >pq2-i k absin(k,c) 1 V a J \ b ) 

1 ” " m, n 

y\ c) • y) cosf^^lsinf^-Wrfy 


. fmtyA 


= ip,? 


(53) 


Since $/-<//?(*', y', c) = /) , expression (53) can be written in component form as 
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= I IV, It 


e 0 £ 0m £ 0» 


Jlyc 

y ^ k, absin(kjC) K ' I ~' V a J"“\ b 

r = rn, n 
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dx'dy' 


(54) 


The total magnetic field inside the enclosure is then obtained from (24) and (54) as 
jUy C _ -j'co v v i v ~ £ o £ ow £ o» f miCyin^f mKX \ rn „( nny } 
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1 1 V,I 


c 0 r = 1 P’ ? 


k t absm(k[C)\ a ) b 
cos (k [Z )I rpqmny 


V 0 J C0S lirJ 


(55) 
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■pqrnny 


(56) 


«f' c = ^IIV,I 


-e 


£ Om £ o« nn, 


m /- = IP’? 

sin (kjz)I r 


, , . ( mnx\ ( nny 

i — t~~Ti — ;ir( _ ^/) cos — cos ~~r 
kj absmikjC) b 1 \ a J \ b 


'pqrnny 


(57) 


For a unique solution the electromagnetic fields in various regions must satisfy continuity condi- 
tions over their common surfaces. The tangential electric fields over the apertures are continuous. 
The tangential magnetic fields over the apertures must also be continuous, thus yielding coupled 


25 



integral equations with the magnetic currents as unknown variables. The coupled integral equa- 
tion in conjunction with the method of moments can be solved for the amplitudes of magnetic cur- 
rents. 


2.3 Derivation of Integral Equation 

The total tangential magnetic fields inside the cavity from apertures on both sides are writ- 


ten as 


H 1 * = + H l ' yi) + H[ Ixc + H'' yc ) 


H 1 ' = (Hy * 0 + H l y Iy ° + Hy XC + H l y yc ) 


Using continuity of tangential magnetic fields across the apertures in the z = 0 plane yields 


H xi 


z = 0 


+h: 


z = 0 


= (H { Ix0 + Hi Iy ° + H Uxc + H Uyc ) 


z = 0 


(58) 


H 


y> 


« +H > 


= (Hy lX0 + Hy y0 + H 1 ’^ + H" yf ) 


z = 0 


(59) 


And matching the tangential magnetic fields across the apertures in z = c plane yields 




= H 


in 


(60) 


(H 1 ^ 0 + Hy y0 + H 1 ^ + H nyC ) 


= H 


III 


(61) 


In deriving equations (58)-(61) it is assumed that the cavity is excited by a plane wave incident 


from z = — °o . 


2.3.1 Integral Equation for Apertures on z = 0 Plane 


Using the field expressions derived in the earlier sections, equation (58) can be written as 
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°> e o U rpq 
2 2 
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(62) 


Now selecting x ¥,.’ p ’ q ’ x as a testing function and use of Galerkin’s method reduces the above equa- 
tion to 
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( 63 ) 


Rearranging the tenns in (63) we get 
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F = V V (1/ Y xlxl +V Y xlyl +4 Y xly2 +R Y xlx2 ' 
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r = 1 P ( i 


( 64 ) 
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Using expression (59) we get 
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Now selecting -O r . , , as a testing function and use of Galerkin’s method yields 
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Rearranging the tenns in (71) we get 
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2.3.2 Integral Equations for Apertures on z = c Plane 
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Now selecting testing function — 'F J .y ljr and use of Galerkin’s method yields 
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ReaiTanging the tenns in (79) we get 
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Rearranging the terms in (84) we get 
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Equations (64), (72), (80), and (87) can be written in a matrix form as 
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The matrix equation (92) can be numerically solved for the unknown amplitudes of equivalent 
magnetic currents induced on the apertures due to given incident field. From the knowledge of 
these amplitudes electromagnetic field inside as well as outside the cavity can be determined. The 
shielding effectiveness of the enclosure with rectangular apertures is then determined from the 
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EM fields using following expression [3] 


E- Shielding(dB) = 
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where Eint is the total electric field at a given point inside the enclosure and Eext is the field at the 
same point in absence of the enclosure. The expressions for electric field components inside the 
cavity due external EM sources are given in Appendix IV. 

3.0 Numerical Results & Discussions 


3.1 Single Aperture Case 

In this section we consider a variety of simple cases and compare the numerical results 
obtained using the present method with results published earlier [3] for validation of the computer 
code. 


Consider a rectangular enclosure with single rectangular aperture on the z = 0 plane and 
illuminated by a plane wave at normal incidence polarized in the x-direction (as shown in Figure 
3). The matrix equation (92) for this case reduces to 
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Using the orthogonality of expansion functions it can be shown = 0, 

hence equation (93) simplifies to 

v xlxl 
_ rpqr'p'q_ 

Furthermore, it has been found from numerical results that the expansion modes p = 1,3, 5,..., 
and q = 0 are strongly excited. Inequation (95) p and q are p = 1, 3, 5 ,... and q = 0 
For a numerical solution of (94) we consider a rectangular cavity with a = 30 cms, 
b = 12 cms, and c = 30 cms with a rectangular aperture of size (10x0.5 ) cms 2 , located at 
( 15, 6 , 0) cms, and illuminated by a plane wave at normal incidence. Assuming only 
p = 1, <7 = 0 expansion mode on the aperture and considering only dominant m = 1, n = 0 
mode inside the cavity, electric field shielding obtained using expression (93) is plotted in Figure 
3.1 along with the results from [3]. It is observed that the numerical data obtained using the 
present method agrees well with the earlier published results. Experimental data from [3] is also 
reproduced in Figure 3.1. 
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Electric Field Shielding (dB) 



Figure 3.1 Electric field shielding at the center of 30 cm x 12 cm x 30 cm enclosure 
wit HO x 0.5 cm 2 aperture located at 15 x 6 cm in z=0 plane with dominant 
cavity mode considered. 


37 




Figure 3.2 Electric field shielding at the center of 30 cm x 12 cm x 30 cm enclosure 
withO x 0.5 cm 2 aperture located at 15 x 6 cm in z=0 plane when 
higher order modes inside cavity and single mode on aperture were considered. 


hr order to study the effect of higher order modes, the electric field shielding is calculated 
using the present method with p = 1 ,q = 0 expansion mode and p' = 1, q' = 0 on the aperture 
and considering m — max = 1 200, n — max = 200 modes including dominant mode inside the 
cavity. The electric field shielding obtained using the present method is plotted in Figure 3.2 
along with the results from [3]. The numerical data in Figure 3.2 shows that consideration of 
higher order modes inside the cavity gives results closer to the measured values. Through numer- 
ical calculations, it was also noted that inclusion of higher order modes on the aperture did not 
alter the electric field shielding values shown in Figure 3.2. 

To validate the present method for a wider aperture, the electric field shielding of a rectan- 
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gular enclosure ( 30 cms x 12 cms x 30 cms) with rectangular aperture (20 cms x 3 cms) as shown 
in Figure 4 is determined as a function of frequency and presented in Figure 4 along with earlier 
published data [3] . It may be observed from the data presented in Figure 4 that numerical results 
obtained by the simple transmission line model used in [3] predicts less electric field shielding 
than the present method. However, the location of minimum shielding predicted by the present 
methods agrees better with the measured data compared to transmission line model of [3]. 


ftfl i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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Figure 4 Electric field shielding at the center of (30 x 12 x 30) cm enclosure with 
(20 x 3.0) cm 2 aperture located at 15 x 6 cm in z=0 plane 


hr the numerical data presented in Figure 4 only a single mode (i.e. p = 1 ,q = 0) over 
the aperture and dominant mode inside the cavity were considered. To study the effects of higher 
order cavity modes on the electric field shielding, a numerical example with the cavity and aper- 
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ture sizes same as shown in Figure 4 is considered. The electric field shielding calculated at the 


center of the rectangular enclosure considering maximum value of m equal to 100 and maximum 
value of n equal to 100 is shown in Figure 5 along with other numerical data. From the numerical 
data presented in Figure 5 it may be concluded that inclusion of higher order cavity modes 
improves electric field shielding estimation. 



Figure 5 Electric field shielding at the center of (30 x 12 x 30) cm" enclosure with 
(20 x 3.0) cm 2 aperture located at 15 x 6 cm in z=0 plane 


Figures 6 and 7 show the electric field shielding when more than a single mode is consid- 
ered at the aperture. The dimensions of the enclosure and the rectangular aperture for Figures 6 
and 7 were the same as shown in Figure 4. From Figures 6 and 7 it may be concluded that for the 
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aperture size and the frequency range considered, a single dominant mode at the aperture is ade- 


quet for achieving convergent results. 



(20 x 3.0) cm 2 aperture located at 15 x 6 cm in z=0 plane 
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Figure 7 Electric field shielding at the center of (30 x 12 x 30) cm 3 enclosure with 
(20 x 3.0) cm 2 aperture located at 15 x 6 cm in z=0 plane 


3.2 Two Aperture Case 

Consider a rectangular enclosure with two rectangular apertures, one on the z = 0 and 
another on the z = c planes. It is assumed that the enclosure is illuminated by a plane wave at 
normal incidence on the aperture at the z = 0 plane. Due to assumed polarization of the inci- 


dence wave it can be shown that V = 0 and A rnn = 0. Hence the matrix equation (92) 
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The matrix equation (96) is solved for U and B using numerical methods. From 
U ...... and B rnn the electric field shielding can be determined using expression (93). For numeri- 

cal solution of (96) we consider a rectangular cavity with a = 30cms, b = 12cms, and 
c = 30 cms with two rectangular apertures. One of the rectangular apertures was of size 
(20x3.0) cm 2 , located at 15, 6, 0 cms, and another rectangular aperture of the same size as the 
first aperture but located at 15, 6, 30 cms. The rectangular cavity is illuminated by plane wave at 
normal incidence to the first aperture. Assuming only p = 1, q = 0 expansion mode on the aper- 
ture and considering only dominant m = 1 , n = 0 mode inside the cavity, the electric field 
shielding obtained using the present method is plotted in Figure 8 (Solid Line). The numerical 
data represented by the solid line in Figure 8 is obtained with the internal coupling between the 

x2x\ xlx2 

apertures not zero (i.e. setting Y r p qr , • • Y “< < < & 0 ). The electric field shielding of the 
rectangular enclosure with two apertures is also obtained by extending the simple transmission 
line model described in [3]. In this model the expression for Z 3 given in equation (6) of reference 
[3] is modified to 

(z ai , + jz g mt K (d-pn> (9r 

3 s Z g + jZ ap tan(k g (d-p)) 

to account for the second aperture. The numerical data obtained using the simple transmission 
line are also plotted in Figure 8. General behaviour of the electric field shielding as a function of 
frequency obtained by the present method and the transmission line method [3] is in good agree- 
ment. It is also observed that when the mutual coupling between the apertures is not considered 
the present method gives shielding results identical to the results obtained using transmission line 
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method. 



Frequency (GHz) 

Figure 8 Electric field shielding at the center of (30 x 12 x 30) cm 3 enclosure with 

two identical apertures with (20 x 3.) cm 2 size. Aperture one located at 
(15 ,6 , 0) cms, aperture two located at (15, 6 , 30) cms. 


In Figure 9 the electric field shielding for the enclosure and aperture dimensions as given 
in Figure 8 is studied when higher order modes inside the cavity and on the apertures are consid- 
ered. 
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Electric Field Shielding (dB ) 


80 



Frequency (GHz) 

Figure 9(a) Effect of higher order aperture inodes on electric field shielding 

at the center of (3() x 12 x 30) cm 3 enclosure withtwo identical apertures 
with (20 x 3.) cnr size. Aperture one located at (15 , 6 , 0) cms, 
aperture two located at (15 , 6 , 30) cms. 
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Frequency (GHz) 

Figure 9(b) Effect of higher order aperture modes on electric field shielding 

at the center of (30 x 12 x 30) cm 3 enclosure withtwo identical apertures 
with (20 x 3.) cur size. Aperture one located at (15 , 6 , 0) cms, 
aperture two located at (15,6, 30) cms. # Cavity modes considered 
(100x100). 


To validate the Modal Method (MM) presented in this report, in the following section, we 
compare the electric field shielding data obtained using the MM method with the FEM-MoM [8] 
techniques. For numerical simulation we consider a rectangular cavity with a = 30 cms, 
b = 1 2 cms, and c = 30 cms with two rectangular apertures. One of the rectangular apertures 
was of size (20x3.0) cm 2 , located at 15, 6, 0 cms, and another rectangular aperture of the same 
size as the first aperture but located at 15, 6, 30 cms. The rectangular cavity is illuminated by a 
plane wave at normal incidence to the first aperture. Assuming two aperture modes and m-max = 
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100, and n-max = 100, the electric field shielding at the center of the cavity is calculated using the 
present method and presented in Figure 10. In Figure 10, the numerical data on electric field 
shielding of the cavity using the hybrid Finite Element Method-Method of Moments (FEM/MoM) 
[ 8] is also presented for various sizes of the apertures. Eexcellent agreement between the results 
obtained from the present and FEM/MoM methods validates the accuracy of the present method. 



Figure 10 Electric field shielding of rectangular enclosure of size (30 x 12 x 30 ) cm3 with 
two identical rectangular apertures ( one on each side) of various sizes when 
illuminated by incident plane wave at normal incidence. 


Figure 1 1 shows a rectangular cavity excited through four apertures; two on each side. The cav- 


ity and aperture sizes used in Figure 1 1 are small enough so that shielding effectiveness of the 


cavity can be calculated using FEM/MoM for validation purpose. The electric field shielding of 


the cavity shown in Figure 1 1 is calculated using the present and FEM/MoM methods and pre- 
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sented in Figure 12. A good agreement between the two results further validate the present 
method. The small disagreement at higher frequencies between the two methods may be attrib- 
uted to the numerical inaccuracy involved in the FEM/MoM methods due to discretization level. 



Figure 1 1 Rectangular cavity with four identical apertures and illuminated by 
plane wave with normal incidence. All dimensions in cm 
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Figure 12 Electric field shielding of rectangular cavity shown in Figure 1 1 . Observation 
point (30,6, 15 ) cm. 

3.3 EM Field Penetration Inside B-757 Aircraft Configuration 

A typical cross section of a B-757 aircraft with various dimensions is shown in Figure 13. 

For simplicity, the inner complexities of passenger seat arrangement and overhead compartments 
are assumed not to be present. The cross section shown in Figure 13 can be approximated by a 
rectangular crosssection as shown in Figure 13. Therefore the entire fuselage of a B-757 can be 
approximated by a rectangular cavity of dimensions shown in Figure 14(a). An external electro- 
magnetic source would cause EM field to penetrate through passenger windows and cockpit win- 
dows into the fuselage. The EM field present inside the fuselage due to any external EM source 
may coupled to sensitive instrumentation wiring and hence cause electronic upset. Furthermore, 
if the EM field strength inside a fuselage is sufficiently large it may cause electric spark on broken 
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wires. It is therefore important to estimate electromagnetic field shielding of a B-757 fuselage or 
for that matter any aircraft passenger cabin. The external EM field may also penetrate a B-757 
fuselage through the cockpit windows. However for an aircraft illuminated from sides passenger 
windows will be dominant contributors. Hence in this report only side windows are considered. 
Approximating the B-757 fuselage by a rectangular cavity with rectangular apertures representing 
passenger windows the computer code developed here is used to determine the field strength 
inside the passenger cabin due to any external EM source. For the worst case scenario, it is 
assumed that the cavity is excited by an electromagnetic plane wave at normal incidence with ver- 
itcal as well as horizontal polarizations. 

For the numerical estimate of electric field shielding of a B-757 fuselage, we first consider 
the effect of two windows, one on each side of a rectangular cavity ( 3600 x 400 x 400 cms) as 
shown in Figure 14(b). The window dimensions selected for the numerical simulation were ( 25 x 

35) cm 2 . The windows were assumed to be located with their centers at (1800 , 200 , 0.0) cms and 
(1800, 200, 400) cms. Using the computer code developed here, the electric field shielding calcu- 
lated at (1800, 120, 200) cms as a function of frequency is presented in Figure 13. In the numer- 
ical simulation it is assumed that the incident wave is polarized in the y-direction. The numbers 

of cavity modes considered in this example were m = 0, 1, 2, ,1200 and 

n = 0,1,2, ....200. Any further addition of cavity modes was found to have no appreciable 
effect on the shielding calculation. For the numerical results shown in Figure 15 only a single 
dominant mode (i.e. p = 1 ,q = 0 was considered on the window apertures. For comparison, 
the measured values of electric field shielding for a B-757 aircraft at two separate frequencies are 
also shown in Figure 15. However, these measurements were done under different incident wave 
conditions and hence are shown here only for illustrating that the numerical data obtained using 
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the present method are in the range of measured values. 
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Figure 13 Geometry of B-757 aircraft configuration with cross sectional and 
longitudinal view. Dotted line shows approximated cross section. 
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Figure 14(a) Closed rectangular cavity approximating B-757 fuselage with 

rectangular apertures on side walls representing passenger windows. 
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400 cms 


Figure 14(b) Closed rectangular cavity approximating B-757 fuselage with 

two rectangular apertures on side walls representing passenger windows. 


Figures 15-17 shows the electric field shielding when higher order modes on the apertures 
were considered, hr Figures 15-17, the incident wave was considered to be parallel to the longer 
dimensions of the aperture windows. For the frequency range shown in Figures 15-17 it is clear 
that single p = 1 ,q = 0 mode on the aperture is adequate for obtaining convergent results. Fig- 
ure 18 shows the electric field shielding of a rectangular cavity with two apertures when the inci- 
dent wave was polarized parallel to the shorter dimension of the window. The dimensions of the 
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cavity and aperture used in numerical data shown in Figure 1 8 were the same as given in Figure 
15. Only a single aperture mode was considered in the numerical data shown in Figure 18. 



Frequency (MHz) 

Figure 15 Electric field shielding of rectangular cavity approximating fuselage of 
B-757 aircraft. Two windows on each side of cavity were considered. 
Cavity size ( 3600 x 400 x 400) cm 3 , window size (25 x 35) cm 2 , window 
locations ( 1800, 200, 0.0) cms and (1800, 200, 400) cms. 

Field point location (1800, 120, 200) cms. 
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Frequency (MHz) 


Figure 16 Electric field shielding of rectangular cavity approximating fuselage of 
B-757 aircraft. Two windows on each side of cavity were considered. 
Dimensions of Cavity, aperture and field point are as given in Figure 15 


55 




Electric Field Shielding (dB ) 
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Figure 18 Electric field shielding (for horizontal polarization) of rectangular cavity 
approximating fuslelage ofB-757 aircraft. Two windows on each side of 
cavity were consideredDimensions of cavity, aperture and field point 
are as given in Figure 15 


3.4 Convergence Test For Large size Cavities 

The numerical estimate of electric field shielding in an electrically large size cavities illu- 
minated through large size the apertures depends upon the number of cavity modes as well as 
aperture modes considered in the numerical simulation. An adequate number of cavity modes and 
aperture modes must be considered to achieve numerical stable results. In this section, depen- 
dence of the electric field shielding on number of cavity modes and aperture modes is studied 
through numerical examples. 

For numerical estimate of electric field shielding, we first consider a rectangular cavity 
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( 5400 x 600 x 600 ) cm , which approximates the fuselage of a Boeing-747. For numerical sim- 
ulation only two windows; one on each side of fuselage were considered. The aperture dimen- 
sions selected for the numerical simulation were ( 25 x 35) cm 2 . The window was assumed to be 
located with their centers at (2700, 300, 0.0) cms and (2700, 300, 600) cms. Using the computer 
code developed here, the electric field in volts/meter calculated at (2700, 300, 300) cms as a func- 
tion of mode index n is presented in Figure 19. In the numerical simulation it is assumed that 
the incident wave is polarized in the y-direction. In Figure 19, stable results are obtained for the 
cavity mode index n higher than 200. For Figure 19, only a single mode at the aperture was con- 
sidered. 

For the cavity and window dimensions as shown in Figure 19, the electric field at the cen- 
ter of cavity is calculated and shown in Figure 20 as a function of cavity mode index m . From 
Figure 20 it is clear that m > 1200 achieves numerically stable results. 

After selecting the cavity mode indices m, n for numerically stable results, the electric 
field at the center of the cavity with dimensions as shown in Figure 19 is calculated as a function 
of aperture modes and shown in Figure 21. From Figure 21 it may be concluded that for normal 
incidence the 4 or 5 aperture modes are sufficient to obtain numerically stable results. Consider- 
ation of higher order modes on the aperture improves the accuracy of the numerical calculation 
but increases the computational time significantly. For fast computation of electric field shielding, 
we consider only single mode at the aperture. 
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Figure 19 Electric field amplitude at (2700, 300, 300) cms inside a rectangular cavity 

( 5400x600x600 ) cm 3 illuminated $y plane wave at normal incidence through 
rectangular windows (25 x 35) crrr as a function of cavity mode index n. 

Other parameters: cavity mode index m = 1200, frequency = 2.87510864 GHz. 
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Figure 20 Electric field amplitude at (2700, 300, 300) cms inside a rectangular cavity 
( 5400x600x600) cm 3 illuminated by plane wave at normal incidence through 
rectangular windows (25 x 35) citT as a function of cavity mode index m. 
Other parameters: cavity mode index n = 200, frequency = 2.87510864 GHz. 
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Figure 21 Electric field amplitude at (2700, 300, 300) cms inside a rectangular cavity 
( 5400x600x600) cm 3 illuminated ^y plane wave at normal incidence through 
rectangular windows (25 x 35) cm“ as a function of number of modes on apertures. 
Other parameters are cavity mode index n = 200, m = 1200, 
frequency = 2. 8745GHz. 
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4.0 Conclusion 


A method suitable for estimation of electric field shielding effectiveness of large but reg- 
ular shaped rectangular cavity with multiple rectangular apertures is presented. Assuming appro- 
priate electric field distribution on the aperture, fields inside the cavity are determined using 
rectangular cavity Green’s function. Electromagnetic fields outside the cavity and scattered due 
to the aperture are obtained using the free space Green’s function. Matching the tangential mag- 
netic field across the apertures, the integral equation with aperture fields as unknown variables is 
obtained. Tire integral equation is solved for unknown aperture fields using the Method of 
Moments. From the aperture fields the electromagnetic shielding effectiveness of the rectangular 
cavity is determined. 

Numerical results on electric field shielding of a rectangular enclosure are validated with 
data available in the literature. Effects of cavity and aperture modes for large size enclosure are 
studied. Approximating a Boeing-757 aircraft cabin by a rectangular cavity and passenger win- 
dows by rectangular aperture electromagnetic field penetrations inside the cabin due to external 
plane wave at normal incidence are estimated. Electromagnetic field penetrations due to horizon- 
tal as well as vertical polarizied incident fields are studied. Though in the present report electro- 
magnetic field penetration through two apertures, one on each side of fuselage is estimated, the 
computer code can be used to simulate effects of multiple apertures on each side of fuselage. 

One advantage of the present method is that the order of matrix equation to be solved in 
the present approach is very small compared to the order of matrix equation one encounters in the 
Finite Element Method. Another advantage is that the computer storge requirement for the 
present approach is negligibly small compared to the Finite Element Method and Finite Differ- 
ence Time Domain method. However, because of the use of the cavity Green’s function it is very 
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difficult to take into account the complex geometries and objects that may be present in the pas- 
senger cabins. Also, the method does not take into account the losses that will be present in prac- 
tical aircraft cabin. Inspite of its inadequateness in handling losses in the system and complex 
geometries, the present method gives quick and fairly accurate estimates of electromagnetic 
shielding effectiveness of rectangular enclosures. 
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In this appendix we evaluate the integral W'VrpqxW’ / ) sin (^)cos {^pjdx'dy' in ; 
closed form. Substituting for ¥ rpqx (x', y') the expression for I rpqmnx can be written as 
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In above equation x r , y r are the coordinates of the center of r aperture. Performing the integra- 
tion expression for I rpqmnx can be written as 
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Appendix II 

hi this appendix we evaluate the integral JJ'Wr.yjcos^jsin 
closed fonn. Substituting for Q> y') the expression for I rpqmny can be written as 
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In above equation x r , y r are the coordinates of the center of r aperture. Perfonrung the integra- 
tion expression for I rpqmny can be written as 
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Appendix III 

The Fourier transforms (|> and \|r of O and *P , respectively are derived in 
this appendix. Using the defination r (k r , k v ) = [['V rnfrv e * kx * ^ kyy dxdy, the Fourier trans- 
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Appendix IV 

In this appendix, the expressions for electric field inside the cavity due to magnetic cur- 
rent sources present on its walls are given. Using the electric vector potential given in (31) and 
expression (23), the electric field components inside the cavity due to the x component of mag- 
netic current sources at z = 0 plane are obtained as 
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Likewise the electric field components due to the x-component of magnetic cuiTent sources at 
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